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PROCEDE D 1 EGALISATION NUMERIQUE, ET RECEPTEUR 
DE RADIOCOMMUNICATION METTANT EN CEUVRE UN TEL PROCEDE 

La presente invention concerne legalisation 
numerique des signaux. Elle trouve une application 
5 importante dans le domaine des radiocommunications . 

Le procede s 1 applique lorsqu' on re<?oit Tin signal 
issu d'un emetteur par 1' intermediaire d'un canal de 
transmission entre emetteur et recepteur, dont la reponse 
est connue oil a ete prealablement estimee . Un probleme 
10 principal qui se pose alors est celui du compromis entre 
les performances de l'egaliseur et sa complex! te. 

Une estimation complete, selon le maximum de 
vrai semblance, de tous les symboles discrets composant le 
signal emis est possible, par exemple en employant 
15 1' algori thine de Viterbi (voir G.D. Forney Jr. : « The 
Viterbi Algorithm », Proc. of the IEEE , Vol. 61, No. 3, 
mars 1973, pages 268-278). Neanmbins, des que la reponse 
impulsionnelle des canaux devient longue ou que le nombre 
de valeurs discretes possibles des symboles devient 
20 important, la complexite exponentielle de ces methodes les 
rend impraticables . 

On considere le cas d'un canal de radiocommunication 
servant a la transmission d'un signal compose de sequences 
ou trames success ives de n symboles d k (l<k<n) . Les 
25 symboles d k sont a valeurs discretes : binaires (±1) dans 
le cas d'une modulation de type BPSK (binary phase shift 
keying), quaternaires (±l±j) dans le cas d'une modulation 
de type QPSK (quaternary phase shift keying) . . . 

Apres conversion en bande de base, numerisation et 
30 filtrage adapte, un vecteur Y du signal re?u refletant les 
symboles emis sur la duree d'une trame a pour expression : 
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ou W+l est la longueur, en nombre de bits, de la repon.se 
impulsionnelle estimee du canal, r == (r Q , t 1 , . . . , r w ) est la 

reponse impulsionnelle estimee du canal, les r q etant des 
5 nombres complexes tels que r^O si q<0 ou q>W, y k est le 
k-ieme echantillon complexe re?u avec l<k<L=n+W, et Y N est 
un vecteur de taille L compose d 1 echantillons de bruit 
additif Y N/k - La reponse impulsionnelle estimee r tient 
compte du canal de propagation, de la mise en forme du 

10 signal par I'emetteur et du filtrage de reception. 

La matrice A de taille Lxn a une structure de type 
Toeplitz le long de sa diagonale principale, c'est-a-dire 
que si a ± ^ designe le terme situe a la i-ieme ligne et a 
la j-ieme colonne de la matrice A, alors <x i+1 ^ j + i=<*i j pour 

15 l<i<L-l et l<j<n-l. Les termes de la matrice A sont donn4s 
par : a x ^=0 pour Kj<n (A n 1 a done que des zeros au-dessus 
de sa diagonale principale) ; a i f i-° pour W+l<i<L 
(structure de matrice -bande) ; et a A ^ r i-i pour l<i<W+l . 

La relation matricielle (1) exprime que le signal 

20 re<?u Y est 1' observation, affectee d'un bruit additif, du 
produit de convolution entre la reponse impulsionnelle du 
canal et les symboles emis . Ce produit de convolution peut 
encore s'exprimer par sa trans forme e en Z : 

Y<Z) = R(Z) .D(Z) + Y N (Z) (2) 

25 ou D(Z), Y(Z), R(Z) et YN(Z) sont les transformers en Z 
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respectives des symboles emis, du signal regru, de la 
reponse impulsionnelle du canal et du bruit : 

□ INCORPORER Equation. 2 □□□ P ( ?^ (3) 

□ INCORPORER Equation. 2 □□□ VfO (4) 
5 □ INCORPORER Equation. 2 □□□ / f ! J (5) 

One solution classique pour resoudre un systeme tel 
que (1) est la methode dite de forgage a zero (« zero 
forcing ») , suivant laquelle on determine le vecteur □ 
INCORPORER Equation. 2 □□□ZF a n composantes continues qui 
10 minimise 1 1 erreur quadratique □ INCORPORER Equation. 2 

Une discretisation des composantes du vecteur □ 
INCORPORER Equation. 2 □□□ZF relative a chaque canal 
intervient ensuite, souvent par le biais d'un decodeur de 

canal. La solution □ INCORPORER Equation. 2 DD ZF au sens 
15 des moindres carres est donnee par : D Z p = (A H A) ~ 1 A H Y , ou A H 
designe la matrice transposee conjuguee de A. On est alors 
raxnene au probleme de I 1 inversion de la matrice 
hermitienne definie positive A H A. Cette inversion peut 
etre realisee par divers algorithmes classiques, d'une 
20 maniere directe (methodes de Gauss, de CholesJcy...) ou par 
des techniques iteratives (algorithmes de Gauss-Seidel , du 
gradient . . . ) . 

L ' erreur d' estimation D— D 2F est egale a (A H A) _1 A H Y N , 
ce qui montre que la solution obtenue est affectee d'un 
25 bruit de variance : 

a 2 = e(||d-D zf || 2 ) = N 0 x Trace[(A H A)" 1 ] (6) 

ou N 0 est la den site spectrale de puissance du bruit. On 
voit qu'il se produit une amplification du bruit (« noise 
enhancement ») quand la matrice A H A est mal condi tionnee , 
30 c'est-a-dire quand elle a une ou plusieurs valeurs propres 
proches de 0 . 

Cette amplification du bruit est le principal 
inconvenient des methodes de resolution classiques. Dans 
la pratique, les cas de mauvais conditionnement de la 



matrice A A sont frequents, particulierement en presence 
de trajets multiples de propagation. 

On connait un moyen relativement simple de remedier 
en par tie a cet inconvenient, en acceptant dans la 
5 solution un residu d' interference, c'est-a-dire en 
adoptant non pas la solution optimale au sens des moindres 
carres, mais la solution: ^mmse~ (A h A+ N q) -1 A h Y , ou N Q 
designe une estimation cLe la densite spectrale du bruit, 
que le recepteur doit alors calculer . Cette methode est 
10 connue sous le nom de MMSE (minimum mean square error) 
Elle permet de diminuer la variante d' estimation par 
rapport a la methode de « zero forcing », mais en 
introduisant un biais . 

Les methodes de « zero forcing » et analogues 
15 reviennent a operer un filtrage inverse du signal re<?u par 
un filtre, modelisant la fonction de transfert 1/R(2) , 
calcule par une certaine approximation (quadratique dans 
le cas du « zero forcing ») . Lorsqu'une ou plusieurs 
racines du polynome R(Z) (equation (5)) sont situees sur 
20 le cercle unite, le filtre inverse theorique presente des 
singularites telles qu' il ne peut pas etre estime par une 
approximation satisf aisante . Dans le cas de 

1 ' approximation quadratique , ceci correspond a la 
divergence de la variance de 1'erreur a 2 lorsque la 
25 matrice A H A a une valeur propre nulle (relation (6)) . 

Ce probleme n' est pas rencontre dans les methodes 
telles que 1' algorithme de Viterbi qui prennent 
intrinsequement en compte la nature discrete des symboles, 
mais qui requierent une puissance de calcul tres 
30 superieure pour les systemes de taille in5>brtante. 

La presente invention a pour but de proposer un 
procede d' egalisation procurant un bon compromis entre la 
fiabilite des estimations et la complexite de l'egaliseur. 

Un autre but est d' obtenir un egaliseur necessitant 
35 une puissance de calcul raisonnable et capable de traiter, 
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avec des performances comparables a celles d'un egaliseur 
de Viterbi, des signaux dont les symboles ont un noxnbre 
d' etats relativement el eves et/ou des signaux regus 
suivant canal de reponse impulsionnelle relativement 

5 etalee. 

L' invention propose ainsi un procede d f egalisation 
nuiuerique , pour es timer des symboles discrets d'un signal 
transmis a partir d 1 echantillons numeriques d'un signal 
regu par 1' intermediaire d'un canal de transmission 
10 represents par une reponse impulsionnelle finie de W+l 
coefficients , W etant un entier plus grand que 1 . Ce 
procede comprend les etapes suivantes : 

- determiner les W racines dans le plan complexe de 
la transformee en Z de la reponse impulsionnelle du 

15 canal ; 

- repartir les W racines en un premier ensemble de 
W-p racines et un second ensemble de p racines , p etant un 
entier plus grand que 0 et plus petit que W, les racines 
du second ensemble etant plus proches du cercle unite que 

20 celles du premier ensemble selon un critere de distance 
determine dans le plan complexe ; 

- obtenir un signal intermediaire en appliquant au 
signal regu une premiere methode d' egalisation sur la base 
d'une reponse impulsionnelle finie dont la transformee en 

25 Z, formee par un polynome en Z" 1 de degre W-p, a pour 
racines les W-p racines du premier ensemble ; et 

- obtenir des estimations des symboles discrets du 
signal transmis en appliquant au signal intermediaire une 
seconde methode d' egalisation sur la base d'une reponse 

30 impulsionnelle finie dont la transformee en Z, formee par 
un polynome en Z" 1 de degre p, a pour racines les p racines 
du second ensemble. 

La « premiere methode d' egalisation » sera 
generalement choisie de fa<?on a traiter les symboles 

35 inconnus comme des variables continues. Elle conduit alors 
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a une operation semblable a un filtrage Inverse dont la 
fonction de transfert serait d' une forme approchant 
1'expression 1/R S (Z) , ou R S (Z) designe le polynome en Z" 1 
de degre W-p ayant pour racines les W-p racines les plus 
5 eloignees du cercle unite . Elle peut notamment etre du 
type « zero forcing ». Cette operation ne genere qu'une 
amplification du bruit reduite, puis que les racines de la 
fonction de transfert en Z associee sont relativexnent 
eloignes du cercle unite. 

1Q Pour les p racines les plus proches du cercle unite, 

on adopte des mesures permettant de s'affranchir ou de 
limiter 1' incidence du probleme de 1' amplification du 
bruit. On peut choisir une methode MMSE ou analogue comxne 
« seconde methode d' egalisation ». Toutefois, cette 

15 seconde methode tiendra avantageusement compte de la 
nature discrete des symboles inconnus . Elle pourra 
notamment reposer sur un algorithme a treillis, tel que 
1' algorithme de Viterbi, dont la mise en ceuvre est 
courante dans les egaliseurs de canal quand la taille du 

20 systeme n'est pas trop grande. 

La seconde methode d' egalisation est generalement 
d'une mise en ceuvre plus complexe que la premiere. Dans 
chaque cas particulier, le choix du nombre p permet de 
rechercher le meilleur compromis entre la fiabilite des 

25 estimations, qui fait preferer les valeurs elevees de p, 
et la complex! te de l'egaliseur^ qui fait preferer les 
valeurs faibles de p. 

Un autre aspect de la presente invention se rapporte 
a un recepteur de radiocommunication r comprenant : 

30 - des moyens de conversion pour produire des 

echantillons numeriques a partir d'un signal radio re<?u 
par 1' intermediaire d'un canal de transmission represents 
par une reponse impulsionnelle finie de W+l coefficients, 
W etant un entier plus grand que 1 ; 

35 - des moyens de mesure de la reponse impulsionnelle 
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du canal ; 

- des moyens de calcul des W racines dans le plan 
complexe de la transformee en Z de la reponse 
impulsionnelle mesuree ; 

5 - des moyens de repartition des W racines en un 

premier ensemble de W-p racines et un second ensemble de p 
racines, p etant un entier plus grand que 0 et plus petit 
que W, les racines du second ensemble etant plus proch.es 
du cercle unite que celles du premier ensemble selon un 
10 critere de distance determine dans le plan complexe ; 

- un premier etage d' egalisation pour obtenir un 
signal intermediaire en appliquant au signal re<?u une 
premiere methode d' egalisation sur la base d' une reponse 
impulsionnelle finie dont la transformee en Z, formee par 

15 un polynome en Z" 1 de degre W-p, a pour racines les W-p 
racines du premier ensemble ; et 

- un second etage d' egalisation pour obtenir des 
estimations de symboles discrets d'un signal transmis sur. 
le canal en appliquant au signal intermediaire une seconde 

20 methode d' egalisation sur la base d'une reponse 
impulsionnelle finie dont la transformee en Z, formee par 
un polynome en Z* 1 de degre p, a pour racines les p racines 
du second ensemble. 

D 1 autres particularites et avantages de la presente 

25 invention apparaitront dans la description ci-apres 
d'exemples de realisation non limitatifs, en reference aux 
dessins annexes, dans lesquels : 

- la figure 1 est un schema synoptique d' un exemple 
de recepteur de radiocommunication selon 1' invention ; 

30 - la figure 2 est un organigramme montrant un mode 

de realisation du procede selon 1' invention ; et 

- la figure 3 est un diagramme illustrant les 
performances du procede. 

Le recepteur represents sur la figure 1 comporte un 
35 etage radio 1 qui regoit le signal radio capte par 
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1' antenne 2 et le convertit en bande de base. Le signal en 
bande de base est numerise par un convertisseur 
analogique-numerique 3, puis fourni a un filtre de 
reception 4. Le filtre 4 assure un filtrage adapt e a la 
5 mise en forme des signaux par l'emetteur. II delivre un 
signal numerique a raison d' un echantillon complexe par 
symbole emis . 

Ce signal numerique est fourni a un demodulateur 
comprenant d'une part un module 6 de synchronisation et 

10 d' estimation de canal, et d' autre part un egaliseur 7. 

La synchronisation et 1' estimation de canal sont par 
exemple effectuees de maniere classique a l'aide d'une 
sequence de synchronisation incluse par l'emetteur dans 
chaque trame de signal. La detection de cette sequence, 

15 connue du recepteur, permet d'une part de synchroniser le 
recepteur par rapport a la structure temporelle des trames 
emises, et d' autre part d'estimer la reponse 
impulsionnelle r = (r Q , r x , . . . , r w ) du canal sur lequel les 

trames sont transmises . La reponse impulsionnelle calculee 
20 par le module 6 est fournie a 1' egaliseur 7. 

L' egaliseur 7 f onctionne par exemple conf ormement a 
1' organigramme represente sur la figure 2 pour traiter 
chaque trame synchronisee du signal re9U, se presentant 



sous la forme d' un vecteur Y= 



avec L=n+W en reprenant 



25 les notations precedentes . 

Le module d' estimation de canal 6 ayant fourni les 
W+l coefficients complexes r q de la reponse impulsionnelle 

estimee du canal, la premiere etape 10 consiste a 
rechercher les W racines de la trans formee en Z de cette 
30 reponse impulsionnelle, donnee par 1' equation (5). 
Diver ses me th odes classiques de recherche de racines 
complexes d' un polynome peuvent etre utilisees a 1' etape 
10. On pourra a cet egard se reporter a l'ouvrage de 



10 



15 



20 



25 
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des Equations 
Type F (x) =0 » , 



E. DURAND : « Solutions Numeriques 
Algebriques ; Tome I : Equations du 
Editions Masson, 1960. 

Les W racines complexes ainsi trouvees a x , a 2 , - . . , 

sont ensuite ordonnees de fa^on a pouvoir les repartir en 
deux ensembles, l'un con tenant les W-p racines les plus 
eloignees du cercle unite, et 1' autre les p racines les 
plus proches du cercle unite . 

Pour cela, une distance 5 q est calculee a l'etape 11 

pour chacune - des racines a q (l<q<W) . Cette distance est 
avantageusement obtenue de la maniere suivante : 



5 q = 



1 - 



1-1/ 



SI 

si 



a, 



a, 



<1 
>1 



(7) 



de la fonction de 



A l'etape 12, les racines <x q 
transfert R(Z) sont triees dans l'ordre des distances 
decroissantes : 5 X > 5 2 ^ ... - &vf- ° n separe alors les W-p 
premieres racines a 1 ,...,a w _ p , qui sont les plus eloignees 
du cercle unite, des p racines restantes oc w _ p+1 , . . . / a w 

A l'etape 13, l'egaliseur 7 developpe un polynome en 
Z" 1 defini par : 

W-p w-p 



R S (Z) 



= n (l - a q . Z- 1 ) = e *< 



(8) 



q=l 



q=0 



Ceci permet de determiner les coefficients s q de la 

fonction de transfert R s (Z) associee a la reponse 
impulsionnelle s = (s 0 , s x , . . . , s w _ p ) d'un canal virtuel, qui 
correspondrait au canal de transmission estime avec 
elimination des contributions les plus proches des zones 
de singularity. 

On peut alors proceder a une premiere egalisation 14 
revenant a effectuer un filtrage inverse approchant la 
fonction de transfert 1/R S (Z). Plusieurs implementations 
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1Q 



15 



peuvent etre retenues pour effectuer ce filtrage inverse. 
On peut notamment effectuer une egalisation par « zero 
forcing » comme indique precedemment . Au sujet de ces 
methodes , on pourra se reporter a l'ouvrage de J.G. 
Proakis : « Digital Communications » McGraw-Hill, 

2 e edition, 1989. 

Le filtrage inverse 14 produit un signal 
intermediaire sous la forme d' un vecteur Y' de L'=n+p 
echantillons y' x , . . . ,y' L , . Dans le cas d'une methode de 

« zero forcing », le vecteur Y' est obtenu par la relation 
matricielle : 

Y' = (A' H A' ) -1 A' H Y (9) 
Dans 1' expression (9), A' designe une matrice de n+W 
lignes et n+p colonnes ayant une structure de Toeplitz, 
formee a partir des coefficients s q du polyndme R S (Z) : 



A' = 



s 0 



»W-p 

0 i 



0 



0 
0 



0 
0 

s 0 



0 
0 



'W-p 

0 s w _ p 



0s w- P ; 



(10) 



20 



Grace au tri des racines <x q , les valeurs propres de 

la matrice A' H A' sont relativement eloignees de 0 . 

En variante, on pourrait realiser le filtrage 
inverse en mettant en cascade W-p cellules de filtrage 
correspondant chacune a 1' inverse d'une fonction de 



transfert r| (Z) = 1 - a q Z -1 , pour l<q<W-p 



Si 



filtre inverse de R q (Z) est irrealisable . Si 
peut developper 1/R q (Z) sous la forme : 



a, 



= 1, le 
< 1 , on 
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— ^— = l + ct CT . Z' 1 + a 2 Z~ 2 + . . . + a* Z" m + . . . (11) 

rJ(z> 

Le developpement (11) est causal, et stable puis que le 
domaine de convergence contient le cercle unite . La 
cellule de filtrage inverse peut done etre realisee sous 
forme transverse ou sous forme recursive. 

> 1, on peut developper 1/r| (Z) sous la 



Si a q 



forme 



* = -a- 1 . Z . (l + a a \ Z + a* Z 2 + . . . + a* Z m + . . .) (12) 
R|(Z) q - V ' 

Ce developpement (12) est anti-causal et stable. Pour la 
realisation de la cellule de filtrage inverse, on tronque 
le developpement (12) et on adopte une implementation sous 
forme transverse. L' anti-causalite provoque un retard 
correspondant a la longueur de la reponse retenue . 

On note que les developpements (11) et (12) 
justifient le critere de distance au cercle unite 8 q 

utilise conformement a la relation (7) . 

A l'etape 15, l'egaliseur 7 developpe un polynome de 
degre p en Z" 1 , dont les racines correspondent aux p 
racines de R(Z) les plus proches du cercle unite, tel que 
R(Z) = R S (Z) .R J (Z) : 

W P 

R^Z) = r 0 . ft (l - a q . Z- 1 ) = I .t q . Z-« (13) 

q=W-p+l q=0 

Les coefficients complexes t q definissent la reponse 
impulsionnelle d'un autre canal de transmission virtuel, 
dont legalisation par une methode de type « zero 
forcing » ou analogue poserait des problemes 
d' amplification du bruit. 

Le signal intermediaire Y' est alors soumis a une 
egalisation selon une autre methode, sur la base de la 
reponse impulsionnelle t = it Q ,t lf . . . ,t p ) . Cette seconde 
egalisation 16 est avantageusement effectuee a l'aide d'un 
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treillis de Viterbi (voir 1' article precite de G.D. Forney 
Jr., ou l'ouvrage precite de J.G. Proakis) . 

Le second etage d' egalisation 16 produit les 
estimations d k des symboles de la trame (l<k<n) . Ces 
estimations d k fonnant la sortie de l'egaliseur 7 peuvent 
®tre fournies a un module de desentrelacement 8 puis a un 
decodeur de canal 9 qui detecte et/ou corrige 
d' eventuelles erreurs de transmission. 

La figure 3 illustre les performances du procede 
dans le cas de la transmission d'une trame de signal selon 
le format du systeme radiotelephonique cellulaire europeen 
GSM, en remplaqant la modulation binaire de type GMSK par 
une modulation de phase a huit etats (modulation 8-PSK) . 
La reponse impulsionnelle du canal etait tronquee a cinq 
15 temps bits (W=4) . La figure 3 montre la dependance entre 
le taux d'erreur binaire BER, exprime en % , et le rapport 
Eb/NO entre l'energie par bit et la densite spectrale du 
bruit, exprime en decibels. Le BER est celui observe dans 
les estimations des symboles apres le desentrelacement et 
20 le decodage de canal effectues conf ormement aux methodes 
employees dans le GSM. La courbe I montre les resultats 
procures par la methode de « zero forcing » pur, c'est-a- 
dire dans le cas limite ou p=0 . La courbe II montre le 
resultat theorique qui serai t obtenu en egalisant le canal 
25 purement avec 1'algorithme de Viterbi (cas limite ou p=W) . 
Dans la pratique, le treillis correspondant devrait 
comporter 8 4 =4096 etats, de sorte que 1' egaliseur de 
Viterbi serait irrealisable avec les techniques actuelles . 
L'ecart entre les courbes I et II illustre la superior! te 
de 1'algorithme de Viterbi qui delivre les estimations 
selon le maximum de vraisemblance . 

Les courbes III et IV montrent les resultats obtenus 
par le procede selon 1' invention, respectivement dans les 
cas ou p=l et p=2 . On voit 1' amelioration tres sensible 
des resultats deja obtenue pour la valeur p=l par rapport 
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au « zero forcing » pur . 

A titre indicatif, 1' egalisation d'tine trame de 
signal GSM par 1' algorithme de Viterbi pur dans les 
conditions de la figure 3 requerrait de l'ordre de 8,45 
5 millions d' operations en virgule flottante, soit environ 
1,83 Gflops, alors que la mise en ceuvre de la presente 
invention dans les memes conditions requiert de l'ordre de 
19000 operations en virgule flottante («4,2 Mflops) dans 
le cas ou p=l, y compris la recherche des racines de R(Z) 
10 et le filtrage inverse 1/R S (Z) par la methode « zero 
forcing ». Ce nombre est de l'ordre de 129000 operations 
(ss28 Mflops) dans le cas ou p=2 , ce qui reste compatible 
avec la puissance des processeurs de traitement de signal 
(DSP) actuellement disponibles . 
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REVENDICATIONS 

1. Procede d' egalisation numerique, pour estimer des 

symboles discrets (d k ) d'un signal transmis a partir 
d'echantillons numeriques (y k ) d'un signal re?u par 
5 l'intermediaire d'un canal de transmission represents par 
une reponse impulsionnelle finie de W+l coefficients 
( r 0' r i' * ■ * ' r w^ ' w etant un entier plus grand que 1, 
comprenant les etapes suivantes : 

- determiner les W racines (a x , <x 2 , . . . f CL^) dans le 

10 plan complexe de la transfonnee en Z (R(2)) de la reponse 
impulsionnelle du canal ; 

- repartir les W racines en un premier ensemble de 
W "P racines (a x , . . . , o^_ p ) et un second ensemble de p 
racines (o^_p +1 / . - - ,o^) , p etant un entier plus grand que 0 

15 et plus petit que W, les racines du second ensemble etant 
plus proches du cercle unite que celles du premier 
ensemble selon un critere de distance determine dans le 
plan complexe ; 

- obtenir un signal intermediaire (Y' ) en appliquant 
20 au signal re?u (Y) une premiere methode d' legalisation sur 

la base d' une reponse impulsionnelle finie dont la 
transfonnee en Z (R s (Z) ) , formee par un polynome en Z" 1 de 
degre W-p, a pour racines les W-p racines du premier 
ensemble ; et 

25 - obtenir des estimations (d^) des symboles discrets 

du signal transmis en appliquant au signal intermediaire 
une seconde methode d' egalisation sur la base d' une 
reponse impulsionnelle finie dont la transfonnee en Z 
(R 1 (Z) ) , formee par un polynome en Z" 1 de degre p, a pour 

30 racines les p racines du second ensemble. 

2 . Procede selon la revendication 1 , dans lequel la 

premiere methode d' egalisation produit le signal 
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intermediaire sous la forme d'un vecteur Y' de n+p 
echantillons (y' 1 , . . . , y' n+p ) obtenu selon la relation : 

Y' = (A /H A' )" 1 A ,H Y 
ou Y est un vecteur forme de n+W echantillons <y k ) du 
signal re<?u, et A' est une ma trice de n+W lignes et n+p 
colonnes ayant une structure de Toeplitz formee a partir 
des coefficients (s q ) dudit polynome en Z _1 de degre W-p 

(R S (Z)) . 

3. Procede selon la revendication 1 ou 2, dans lequel 

la seconde methode d' egalisation comporte la raise en ceuvre 
d'un algorithme de Viterbi. 



4. Procede selon 1'une quelconque des revendications 

1 a 3, dans lequel le critere de distance au cercle unite, 
utilise pour repartir les W racines a^...,^ de la 
15 transformee en Z (R(Z) ) de la reponse impulsionnelle du 
canal entre les premier et second ensembles, s'exprime par 



si ot q 



une distance 5 q de la forme 8 q - 1 — |a q 
de la forme 8 q = 1 - l/|<*q| si |cc q | > 1, pour l<q<W 



< 1 , et 



5. Recepteur de radiocommunication, comprenant : 

20 - des moyens de conversion (1,3,4) pour produire des 

echantillons numeriques (y k ) a partir d'un signal radio 
re<?u par 1' intermediaire d'un canal de transmission 
represents par une reponse impulsionnelle finie de W+l 
coefficients (r 0 , r x , . . . , r w ) , W etant un entier plus grand 

25 que 1 ; 

- des moyens (6) de mesure de la reponse 
impulsionnelle du canal ; 

- des moyens de calcul des W racines (o^ , a 2 , . - . / ct^) 
dans le plan complexe de la transformee en Z (R(Z)) de la 

30 reponse impulsionnelle mesuree ; 




- 16 - 

- des moyens de repartition des W racines en un 
premier ensemble de W-p racines (a x , . . . , a^p) et un second 
ensemble de p racines (o^ +1 , . . . ,0^) , p etant un entier 
plus grand que 0 et plus petit que W, les racines du 

5 second ensemble etant plus proches du cercle unite que 
celles du premier ensemble selon un critere de distance 
determine dans le plan complexe ; 

- un premier etage d' egalisation pour obtenir un 
signal intermediaire en appliquant au signal re<?u (y k ) une 

10 premiere method© d' egalisation sur la base d' une reponse 
impulsionnelle finie dont la transfozmee en 2 <R S (Z)), 
formee par un polynome en Z" 1 de degre W-p, a pour racines 
les W-p racines du premier ensemble ; et 

- un second etage d' egalisation pour obtenir des 
15 estimations ( d k ) de symboles discrets d' un signal transmis 

sur le canal en appliquant au signal intermediaire une 
seconde methode d' egalisation sur la base d' une reponse 
impulsionnelle finie dont la transformee en Z (R 1 ^)), 
formee par un polynome en Z" 1 de degre p, a pour racines 
20 les p racines du second ensemble . 

6. Recepteur selon la revendication 5, dans lequel le 

premier etage d' egalisation est agence pour produire le 
signal intermediaire sous la forme d' un vecteur Y' de n+p 
echantillons (Y' i / • • • / Y' n+p ) obtenu selon la relation : 

25 Y' = (A' H A' )" 1 A' H Y 

ou Y est un vecteur forme de n+W echantillons (y k ) du 
signal regu, et A' est une matrice de n+W lignes et n+p 
colonnes ayant une structure de Toeplitz formee a partir 
des coefficients (s q ) dudit polynome en Z" 1 de degre W-p 

30 (R S (Z)). 
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7. Recepteur selon la revendication 5 ou 6, dans 

lec[uel le second etage d' egalisation est agence pour 
mettre en oeuvre un algorithme de Viterbi. 
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8. Recepteur selon 1'une quelconque des 

revendications 5 a 7, dans lequel les moyens de 
repartition des racines utilisent, pour repartir les W 
racines entre les premier et second ensembles , un critere 
de distance au cercle unite s'exprimant par une distance 5 



de la forme 

5 CT = 1 - l/\a a 



5 q = 1- 



sx 



a, 



< 1 , et de la forme 



sx 



> 1 , pour l<q<W . 
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presence invention concerne 1'egalisation 

• ~ H« a cionaux Elle trouve une application 
numerique des signaux . 

importante dans le doaaine des radiocommunications . 

Le precede s ' applique lorsqu' on recoit un signal 
issu d'un emetteur par 1' intermediate d'un canal de 
transmission entre emetteur et recepteur, dont la reponse 
est connue ou a ete prealablement estimee. Un probleme 
principal qai se pose alors est celui du compromis entre 
les performances de l'egaliseur et sa complexite. 

Une estimation complete, selon le maximum de 
vraisemblance, de tous les symboles discrets composant le 
signal emis est possible, par example en employant 
1-algorithme de Viterbi (voir G.D. Forney Jr. : « The 
Viterbi Algorithm », Proc. of the IEEE , Vol. 61, No. 3, 
mars 1973, pages 268-278). Neanmoins, des que la reponse 
impulsionnelle des canaux devient longue ou que le nombre 
de valeurs discretes possibles des symboles devxent 
important, la complexite exponentielle de ces methodes les 

rend impraticables . 

On considere le cas d'un canal de radiocommunication 
servant a la transmission d'un signal compose de sequences 
ou trames successives de n symboles d k (l<k*n) . Les 
symboles d* sont a valeurs discretes : binaires (±1) dans 
le cas d'une modulation de type BPSK (binary phase shift 
keying), quaternaires (±l±j) dans le cas d'une modulation 
de type QPSK (quaternary phase shift keying) . . . 

Apres conversion en bande de base, numerisation et 
filtrage adapte, un vecteur Y du signal recu refletant les 
symboles emis sur la duree d'une trame a pour expression : 
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r x r 0 0 
r l 

r w 

0 r w 
0 0 r w 



0 
0 

*1 



0 0 r w J 



d 2 



Yn,2 



Yn,Jc 



= A.D + Yi 



N 



(1) 



ou W+l est: la longueur, en n ombre de bits, de la reponse 



impulsionnelle estimee du canal, r = (r 



,r w ) est la 



reponse impulsionnelle estimee du canal, les r q etant des 

5 nombres complexes tels que r q =0 si q<0 ou q>W, y k est le 

k-ieme echantillon complexe regu avec l<k<L=n+W, et Y N est 

un vecteur de taille L compose d 1 echantillons de bruit 
additif y N k - La reponse impulsionnelle estimee r tient 

compte du canal de propagation, de la mise en forme du 
10 signal par l'emetteur et du filtrage de reception. 

La , ma trice A de taille Lxn a une structure de type 
Toeplitz le long de sa diagonale principale, c'est-a-dire 
que si <*± f j designe le terme situe a la i-ieme ligne et a 
la j-ieme colonne de la ma trice A, alors ot i+1 ^ j+i—o^ j pour 

15 l<i<L-l et l<j<n-l . Les tezmes de la matrice A sont donnes 
par : a x j=0 pour Kj<n (A n'a done que des zeros au-dessus 

de sa diagonale principale) ; a i,i = 0 pour W+Ki<L 
(structure de matrice -bande) ; et x™ x L~X P our lSi£W+l . 

La relation matricielle (1) exprime que le signal 
20 re^u Y est 1' observation, af fee tee d' un bruit additif, du 
produit de convolution entre la reponse impulsionnelle du 
canal et les symboles emis . Ce produit de convolution peut 
encore s 7 exprime r par sa trans formee en Z : 

Y(Z) = R(Z).D(Z) + Y N (Z) (2) 

25 ou D(Z), Y(Z), R(Z) et Y N (Z) sont les transformees en Z 
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respectives des symboles emis, du signal re<?u, de la 
reponse ixnpulsionnelle du canal et du bruit : 

n 

D(Z) = 2 d k .Z" k (3) 
Jc=l 

L 

Y(Z) = 2 Yk- Z" k < 4) 
5 R(Z) = ir q .Z^ (5) 

q=0 

Une solution classique pour resoudre un systeme tel 
que (1) est la methode dite de for<?age a zero (« zero 
forcing ») , suivant laquelle on determine le vecteur D 2F a 
n composantes continues qui minimise 1 1 erreur quadra tique 
10 e = ||aD-y|| 2 . Une discretisation des composantes du vecteur 
D relative a chaque canal intervient ensui te 9 souvent 

ZF 

par le biais d'un decodeur de canal. La solution D ZF au 

sens des moindres carres est donnee par : f> Z p = (A H A) " 1 A H Y f 

ou A H designe la matrice transposee conjuguee de A. On est 
15 alors ramene au probleme de 1 1 inversion de la matrice 
hermitienne definie positive A H A. Cette inversion peut 
etre realisee par divers algorithmes classiques, d'une 
maniere directe (methodes de Gauss, de Cholesky...) ou par 
des techniques iteratives (algorithmes de Gauss-Seidel , - du 
20 gradient. . .) . 

L' erreur d 1 estimation D-D 2F est egale a t**K) -laPr^ , 
ce qui montre que la solution obtenue est affectee d'un 
bruit de variance : 

a 2 = e(||d-D zf | 2 ) = N 0 x Trace[(A H A)~ X ] (6) 
25 ou N 0 est la densite spectrale de puissance du bruit. On 
voit qu'il se produit une amplification du bruit (« noise 
enhancement ») quand la matrice A H A est mal conditionnee , 
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c' est-a-dire quand elle a une ou plusieurs valeurs propres 
proches de 0 , 

Cette amplification du bruit est le principal 
inconvenient des methodes de resolution classigues. Dans 
5 la pratique, les cas de mauvais conditionnement de la 
ma trice A H A sont frequents, particulierement en presence 
de trajets multiples de propagation. 

On connait un moyen relativement simple de remedier 
en partie a cet inconvenient, en acceptant dans la 
10 solution un residu d' interference, c'est-a-dire en 
adoptant non pas la solution optimale au sens des moindres 
carres, mais la solution: ^mmse = (^ h A+ N 0 ) -1 A h Y , ou N Q 
designe une estimation de la densite spectrale du bruit, 
que le recepteur doit alors calculer. Cette methode est 
15 connue sous le nom de MMSE (minimum mean square error) 
Elle permet de diminuer la variante d' estimation par 
rapport a la methode de « zero forcing », mais en 
introduisant un biais . 

Les methodes de « zero forcing » et analogues 
reviennent a operer un filtrage inverse du signal regu par 
un filtre, modelisant la fonction de transfert 1/R(Z) , 
calcule par une certaine approximation (quadratique dans 
le cas du « zero forcing ») . Lorsqu'une ou plusieurs 
racines du polynome R(Z) (equation (5)) sont situees sur 
le cercle unite, le filtre inverse theorique presente des 
singularites telles qu'il ne peut pas etre estime par une 
approximation satisf aisante . Dans le cas de 

1' approximation quadratique, ceci correspond a la 
divergence de la variance de I'erreur a 2 lorsque la 
30 ma trice A H A a une valeur propre nulle (relation (6)) . 

Ce probleme n' est pas rencontre dans les methodes 
telles que 1'algorithme de Viterbi qui prennent 
intrinsequement en compte la nature discrete des symboles, 
mais qui requierent une puissance de calcul tres 
35 superieure pour les systemes de taille importante. 
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La presente invention a pour but de proposer un 
procede d' egalisation procurant un bon compromis entre la 
fiabilite des estimations et la complexite de l'egaliseur. 

Un autre but est d' obtenir un egaliseur necessitant 
5 une puissance de calcul raisonnable et capable de traiter, 
avec des performances comparable s a celles d'un egaliseur 
de Viterbi, des signaux dont les symboles ont un nombre 
d' etats relativement eleves et/ou des signaux reqrus 
suivant un canal de reponse impulsionnelle relativement 
10 etalee. 

I*' invention propose ainsi un procede d f egalisation 
numerique, pour es timer des symboles discrets d'un signal 
transmis a partir d 1 echantillons numeriques d'un signal 
re?u par 1' intermediaire d'un canal de transmission 
15 represents par une reponse impulsionnelle finie de W+l 
coefficients , W etant un en tier plus grand que 1 - Ce 
procede comprend les e tapes suivantes : 

- determiner les W racines dans le plan complexe de 
la transformee en Z de la reponse impulsionnelle du 

20 canal ; 

- repartir les W racines en un premier ensemble de 
W-p racines et un second ensemble de p racines, p etant un 
en tier plus grand que 0 et plus petit que W, les racines 
du second ensemble etant plus proches du cercle unite que 

25 celles du premier ensemble selon un critere de distance 
determine dans le plan complexe ; 

- obtenir un signal intermediaire en appliquant au 
signal re£u une premiere methode d' egalisation sur la base 
d' une reponse impulsionnelle finie dont la -transformee en 

30 Z, formee par un polynome en Z" 1 de degre W-p, a pour 
racines les W-p racines du premier ensemble ; et 

- obtenir des estimations des symboles discrets du 
signal transmis en appliquant au signal intermediaire une 
seconde methode d' egalisation sur la base d'une reponse 

35 impulsionnelle finie dont la transformee en Z, formee par 
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un polynome en Z" 1 de degre p, a pour racines les p racines 
du second ensemble . 

La « premiere methode d' egalisation » sera 
generalement choisie de fa<?on a traiter les symboles 
5 inconnus come des variables continues . Elle conduit alors 
a une operation semblable a un filtrage inverse dont la 
f onction de transf ert serait d' une forme approchant 
1' expression 1/R S (Z), ou R S (Z) designe le polynome en Z" 1 
de degre W-p ayant pour racines les W-p racines les plus 
10 eloignees du cercle unite. Elle peut notaxnment etre du 
type « zero forcing ». Cette operation ne gene re qu' une 
amplification du bruit reduite, puis que les racines de la 
f onction de transf ert en Z associee sont relativement 
eloignees du cercle unite. 
15 Pour les p racines les plus proches du cercle unite, 

on adopte des mesures permettant de s'affranchir ou de 
limiter 1' incidence du probleme de 1' amplification du 
bruit. On peut choisir une methode MMSE ou analogue coznme 
« seconde methode d' egalisation ». Toutefois, cette 
20 seconde methode tiendra avantageusement compte de la 
nature discrete des symboles inconnus. Elle pourra 
notamment reposer sur un algorithme a treillis , tel que 
1' algorithme de Viterbi, dont la mise en oeuvre est 
courante dans les egaliseurs de canal guand la taille du 
25 systeme n' est pas trop grande. 

La seconde methode d' egalisation est generalement 
d'une mise en oeuvre plus complexe que la premiere. Dans 
chaque cas particulier, le choix du noinbre p permet de 
rechercher le meilleur compromis entre la fiabilite des 
30 estimations, qui fait preferer les valeurs elevees de p, 
et la complexite de 1'egaliseur, qui fait preferer les 
valeurs faibles de p. 

Un autre aspect de la presente invention se rapporte 
a un recepteur de radiocoznmunication , comprenant : 
35 " des moyens de conversion pour produire des 
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echantillons numeriques a partir d'un signal radio re<?u 
par 1' intenaediaire d' un canal de transmission represents 
par une reponse impulsionnelle finie de W+l coefficients, 
W etant un entier plus grand que 1 ; 
5 - des moyens de mesure de la reponse impulsionnelle 

du canal ; 

- des moyens de calcul des W racines dans le plan 
complexe de la transformee en Z de la reponse 
impulsionnelle mesuree ; 

!0 - des moyens de repartition des W racines en un 

premier ensemble de W-p racines et un second ensemble de p 
racines, p etant un entier plus grand que 0 et plus petit 
que W, les racines du second ensemble etant plus proches 
du cercle unite que celles du premier ensemble selon un 

15 critere de distance determine dans le plan complexe ; 

- un premier etage d' egalisation pour obtenir un 
signal intermediaire en appliquant au signal re<?u une 
premiere methode d' egalisation sur la base d' une reponse 
impulsionnelle finie dont la transformee en Z, formee par 

20 un polynome en Z" 1 de degre W-p, a pour racines les W-p 
racines du premier ensemble ; et 

- un second etage d' egalisation pour obtenir des 
estimations de symboles discrets d'un signal transmis sur 
le canal en appliquant au signal intermediaire une seconde 

_ 25 methode d' egalisation sur -la—base d!_une reponse... 

impulsionnelle finie dont la transformee en Z, formee par 
un polynome en Z' 1 de degre p, a pour racines les p racines 
du second ensemble. 

D'autres particular! tes et avantages de la presente 
30 invention apparaitront dans la description ci-apres 
d'exemples de realisation non limitatifs, en reference aux 
dessins annexes, dans lesquels : 

- la figure 1 est un schema synoptique d'un exemple 
de recepteur de radiocommunication selon 1' invention ; 

35 - la figure 2 est un organigramme montrant un mode 



- 8 - 

de realisation du procede selon 1' invention ; et 

~ la figure 3 est un diagramme illustrant les 
performances du procede . 

Le recepteur represente sur la figure 1 comporte un 
5 etage radio 1 qui re^oit le signal radio capte par 
l'antenne 2 et le convertit en bande de base. Le signal en 
bande de base est numerise par un convertisseur 
analogique-numerique 3, puis fourni a un filtre de 
reception 4. Le filtre 4 assure un filtrage adapte a la 
10 mise en forme des signaux par l'emetteur. II delivre un 
signal numerique a raison d' un echantillon complexe par 
symbol e emis . 

Ce signal numerique est fourni a un demodulateur 
comprenant d' une part un module 6 de synchronisation et 

15 d' estimation de canal, et d' autre part un egaliseur 7. 

La synchronisation et 1' estimation de canal sont par 
exemple effectuees de maniere classique a l'aide d'une 
sequence de synchronisation incluse par l'emetteur dans 
chaque trame de signal. La detection de cette sequence, 

20 connue du recepteur, permet d' une part de synchroniser le 
recepteur par rapport a la structure temporelle des trames 
emises, et d' autre part d' es timer la reponse 
impulsionnelle r = (r Q/ r x , . . . , r w ) du canal sur lequel les 
trames sont transmises. La reponse impulsionnelle calculee 

25 par le module 6 est fournie a 1' egaliseur 1. 

L' egaliseur 7 fonctionne par exemple conformement a 
l'organigramme represente sur la figure 2 pour traiter 
chaque trame synchronisee du signal re?u f se presentant 



sous la forme d' un vecteur Y= 



avec L=n+W en reprenant 



30 les notations precedentes . 

Le module d' estimation de canal 6 ayant fourni les 
W+l coefficients complexes r q de la reponse impulsionnelle 
estimee du canal, la premiere etape 10 consiste a 
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rechercher les W racines de la transformee en Z de cette 
reponse impulsionnelle, donnee par 1' equation (5). 
Diverses methodes classiques de recherche de racines 
complexes d' un polynome peuvent etre utilisees a 1'etape 
5 10. On pour r a a cet egard se reporter a l'ouvrage de 
E . DURAND : « Solutions Numeriques des Equations 
Algebriques ; Tome I : Equations du Type F(x) =0 » , 
Editions Massozi, 1960. 

Les W racines complexes ainsi trouvees a 1 ,a 2/ . . . /Oc^ 

10 sont ensuite ordonnees de fa?on a pouvoir les repartir en 
deux ensembles , l'un con tenant les W-p racines les plus 
eloignees du cercle unite, et 1' autre les p racines les 
plus proches du cercle unite . 

Pour cela, une distance S q est calculee a 1'etape 11 

15 pour chacune des racines a q (l<q<W) . Cette distance est 
avantageusement obtenue de la maniere suivante : 

1 - a J si 
q [l - l/|a q | si 

A 1'etape 12, les racines <x q de la fonction de 
transfert R(Z) sont triees dans l'ordre des distances 
20 decroissantes : 8 1 > 8 2 ^ . - • - $w- ° n separe alors les W-p 
premieres racines a x , . . . ,ct w- p, qui sont les plus eloignees 
du cercle unite, des p" racines res tan tes a w-p+1 , ."T-ya^. 

A 1'etape 13, 1' egaliseur 7 developpe un polynome en 
Z" 1 defini par : 

W-p W-p 

25 R S (Z) = 0 ( X " a q- Z" 1 ) = Z S q- Z " q < 8 > 

q=l q=0 

Ceci permet de determiner les coefficients s q de la 

fonction de transfert R S (Z) associee a la reponse 
impulsionnelle s = (s Q , s x , . . • , s w _ p ) d'un canal virtuel, qui 

correspondrait au canal de transmission estime avec 



a, 



a, 



<1 
>1 



(7) 
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elimination des contributions les plus proches des zones 
de singularity . 

On peut alors proceder a une premiere egalisation 14 
revenant a effectuer un filtrage inverse approchant la 
fonction de transfert 1/R S (Z) . Plusieurs implementations 
peuvent etre retenues pour effectuer ce filtrage inverse. 
On peut notamment effectuer une egalisation par « zero 
forcing » comme indique precedemment . Au sujet de ces 
methodes, on pourra se reporter a I'ouvrage de J.G. 
Proakis : « Digital Communications » McGraw-Hill, 

2 e edition, 1989. 

Le filtrage inverse 14 produit un signal 
intermediaire sous la forme d'un vecteur Y' de L'=n+p 
echantillons y' x , . . . ,y' L , . Dans le cas d' une methode de 

« zero forcing », le vecteur Y' est obtenu par la relation 
matricielle : 



Y' = (A /H A' )~ ± A'* 1 Y (9> 
Dans 1' expression (9), A' designe une matrice de n+W 
lignes et n+p colonnes ayant une structure de Toeplitz, 
formee a partir des coefficients s du polynome R S (Z) : 



-l 



H 



A' = 



s 0 



0 

s 0 

s l 



0 
0 

s 0 

s l 



0 
0 



s W-p 

0 sw-i 



0 0 
Grace au tri des racines a 



0 
0 

s o 



(10) 



* q , les valeurs propres de 
la matrice A' H A' sont relativement eloignees de 0 . 

En variante, on pourrait realiser le filtrage 
inverse en mettant en cascade W-p cellules de filtrage 
correspondant chacune a 1' inverse d'une fonction de 
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<x q | = 1 , le 
aJ < 1, on 



transfert r| (Z) = 1 - a q Z" x , pour l<q<W-p . Si 

filtre Inverse de R q (Z) est irrealisable. Si 

peut developper 1/R q (Z) sous la forme : 

* = 1 + a a . Z" 1 + a* Z' 2 + . . . + a* 2" m + . . . (11) 
rJ(Z) q q q 

Le developpement (11) est causal, et stable puisque le 
domaine de convergence contient le cercle unite . La 
cellule de filtrage inverse peut done etre realisee sous 
forme transverse ou sous forme recursive . 



Si 



a q 



> 1, on peut developper 1/R-(Z) sous la 



10 forme 



R|(Z) q 



1 -1 / —1 —2 5 — m m \ 

= -a a . Z . Il + a q . Z + a q - 2r + . . . + a q . Z m + . . .1 (12) 



Ce developpement (12) est anti-causal et stable. Pour la 
realisation de la cellule de filtrage inverse, on tronque 
le developpement (12) et on adopte une implementation sous 
15 forme transverse. L' anti-causalite provoque un retard 
correspondent a la longueur de la reponse retenue. 

On note que les developpements (11) et (12) 
justifient le critere de distance au cercle unite 8 q 

utilise conf ormement a la relation (7) . 
20 A l'etape 15, l'egaliseur 7 developpe un polynome de 

degre p en Z" 1 , dont les racines correspondent aux p 
racines de R(Z) les plus proches du cercle unite, tel que 
R(Z) = R S (Z) .R 1 (Z) : 

W p 

R J (Z) = r 0 . H(l - <x q . Z' 1 ) = X t q . Z~* (13) 

q=W-p+l q=0 

25 Les coefficients complexes t q definissent la reponse 

impulsionnelle d' un autre canal de transmission virtuel, 
dont 1 ' egalisation par une methode de type « zero 
forcing » ou analogue poserait des problemes 
d' amplification du bruit. 
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Le signal intermediaire Y' est alors soumis a une 
egalisation selon une autre methode, sur la base de la 
reponse impulsionnelle t = (t Q , t x , . . . , t p ) . Cette seconde 
egalisation 16 est avantageusement effectuee a 1'aide d'un 
5 treillis de Viterbi (voir 1' article precite de G.D . Forney 
Jr., ou l'ouvrage precite de J.G. Proakis) . 

Le second etage d' egalisation 16 produit les 
estimations d k des symboles de la trame (l<k<n) . Ces 
estimations d k formant la sortie de 1'egaliseur 7 peuvent 
10 etre fournies a un module de desentrelacement 8 puis a un 
decodeur de canal 9 qui detecte et/ou corrige 
d' eventuelles erreurs de transmission. 

La figure 3 illustre les performances du procede 
dans le cas de la transmission d' une trame de signal selon 
!5 le format du systeme radiotelephonique cellulaire europeen 
GSM, en remplagant la modulation binaire de type GMSK par 
une modulation de phase a huit etats (modulation 8-PSK) . 
La reponse impulsionnelle du canal etait tronqpiee a cinq 
temps bits (W=4) . La figure 3 montre la dependance entre 
20 le taux d'erreur binaire BER, exprime en %, et le rapport 
Eb/NO entre 1' energie par bit et la densite spectrale du 
fc^ruit, exprime en decibels. Le BER est celui observe dans 
les estimations des symboles apres le desentrelacement et 
le decodage de canal effectues conformement aux methodes 
25 employees dans le GSM. La courbe I montre les resultats 
procures par la methode de « zero forcing » pur, c'est-a- 
dire dans le cas limite ou p=0 . La courbe II montre le 
resultat theorique qui serait obtenu en egalisant le canal 
purement avec l'algorithme de Viterbi (cas limite ou p=W) . 
30 Dans la pratique, le treillis correspondant devrait 
comporter 8 4 =4096 etats, de sorte que 1'egaliseur de 
Viterbi serait irrealisable avec les techniques actuelles . 
L'ecart entre les courbes I et II illustre la superiority 
de l'algorithme de Viterbi qui delivre les estimations 
35 selon le maximum de vrai semblance . 
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Les courbes III et IV montrent les resultats obtenus 
par le procede selon 1' invention, respectivement dans les 
cas ou p=l et p=2 . On voit 1' amelioration tres sensible 
des resultats deja obtenue pour la valeur p=l par rapport 
5 au « zero forcing » pur. 

A titre indicatif, legalisation d'une trame de 
signal GSM par 1' algorithme de Viterbi pur dans les 
conditions de la figure 3 requerrait de l'ordre de 8,45 
millions d' operations en virgule flottante, soit environ 

10 1,83 Gflops, alors que la mise en ceuvre de la presente 
invention dans les memes conditions requiert de l'ordre de 
19000 operations en virgule flottante («4,2 Mflops) dans 
le cas ou p=l, y compris la recherche des racines de R(Z) 
et le f iltrage inverse 1/R S (Z) par la methode « zero 

15 forcing ». Ce nombre est de l'ordre de 129000 operations 
(«28 Mflops) dans le cas ou p=2 , ce qui reste compatible 
avec la puissance des processeurs de traitement de signal 
(DSP) actuellement disponibles . 
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REVENDICAT IONS 

1. Procede d 1 egalisation mimerique , pour estimer des 

symboles discrets (d k ) d'un signal transmis a partir 
d' echantillons numeriques (Y^) d'un signal regu par 
5 1' intermediaire d' un canal de transmission represente par 
une reponse impulsionnelle finie de W+l coefficients 
(r 0 ,r x , . . . ,r w ) , W etant un en tier plus grand que 1, 
comprenant les e tapes sulvantes : 

- determiner les W racines (a 1 ,a 2 , . . • ,0^) dans le 

10 plan complexe de la transf ormee en Z (R(Z)) de la reponse 
Impulsionnelle du canal ; 

- repartir les W racines en un premier ensemble de 
W-p racines (a x , . . . , cx^^) et un second ensemble de p 

racines (a^p^ , . . . , o^) , p etant un entier plus grand que 0 

15 et plus petit que W, les racines du second ensemble etant 
plus proches du cercle unite que celles du premier 
ensemble selon un critere de distance determine dans le 
plan complexe ; 

- obtenir un signal intermediaire (Y' ) en appliquant 
20 au signal regu (Y) une premiere methode d' egalisation sur 

la base d' une reponse impulsionnelle finie dont la 
transf ormee en Z (R s (Z) ) , f ormee par un polynome en Z" 1 de 
degre W-p, a pour racines les W-p racines du premier 
ensemble ; et 

25 - obtenir des estimations (d k ) des symboles discrets 

du signal transmis en appliquant au signal intermediaire 
une seconde methode d' egalisation sur la base d'une 
reponse impulsionnelle finie dont la transformee en Z 
(R Z (Z)), f ormee par un polynome en Z' 1 de degre p, a pour 

30 racines les p racines du second ensemble. 

2. Procede selon la revendication 1, dans lequel la 

premiere methode d' egalisation produit le signal 
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intermediaire sous la forme d'un vecteur Y' de n+p 
echantillons (y' 1 , . . . ,y' n+p ) obtenu selon la relation : 

Y' = <A' H A')" 1 A' H Y 
ou Y est un vecteur forme de n+W echantillons (y k ) du 

5 signal re^u, et A' est une ma trice de n+W lignes et n+p 
colonnes ayant une structure de Toeplitz formee a partir 
des coefficients (s q ) dudit polyndme en Z" 1 de degre W-p 

<R S (Z)). 

3 . Procede selon la revendication 1 ou 2 , dans lequel 
10 la seconde methode d' egalisation comporte la mise en oeuvre 

d' un algorithme de Viterbi . 

4. Procede selon l'une quelcongue des revendications 
1 a 3, dans lequel le critere de distance au cercle unite, 
utilise pour repartir les W racines <x lf . . . ,0^ de la 

15 transformee en Z (R(Z)) de la response impulsionnelle du 
canal entre les premier et second ensembles, s'exprime par 

< 1, et 



une distance 5 q de la forme 6 q = 1 — 



a q 



SI 



a q 



de la forme 8 q =1-1/ 



a q 



sx 



a- > 1 , pour l<q<W. 



5. Recepteur de radiocommunication, comprenant : 

20 - des moyens de conversion (1,3,4) pour produire des 

echantillons numeriques (Y^) partir d'un signal radio 

re?u par 1' intermediaire d' un canal de transmission 
represents par une reponse impulsionnelle finie de W+l 
coefficients (r Q , r ± , . . . , r w ) , W etant un entier plus grand 

25 que 1 ; 

- des moyens (6) de mesure de la reponse 
impulsionnelle du canal ; 

- des moyens de calcul des W racines (a x , a 2 , - . . , o^) 

dans le plan complexe de la transformee en Z (R(Z)) de la 
30 reponse impulsionnelle mesuree ; 
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- des moyens de repartition des W racines en un 
premier ensemble de W-p racines (a x , . . . , cc^p) et un second 
ensemble de p racines (a^^ , . . . , o^) , p etant un entier 
plus grand que 0 et plus petit que W, les racines du 

5 second ensemble etant plus proches du cercle unite que 
celles du premier ensemble selon un critere de distance 
determine dans le plan complexe ; 

- un premier etage d' egalisation pour obtenir un 
signal intermediaire en appliquant au signal regu (y k ) une 

10 premiere methode d' egalisation sur la base d' une reponse 
impulsionnelle finie dont la transformee en Z (R S (Z)), 
formee par un polynome en Z" 1 de degre W-p, a pour racines 
les W-p racines du premier ensemble ; et 

- un second etage d' egalisation pour obtenir des 
15 estimations (d k ) de symboles discrets d' un signal transmis 

sur le canal en appliquant au signal intermediaire une 
seconde methode d' egalisation sur la base d' une reponse 
impulsionnelle finie dont la transformee en Z (R 1 ^)), 
formee par un polynome en Z" 1 de degre p, a pour racines 
20 les p racines du second ensemble. 

6 . Recepteur selon la revendication 5 , dans lequel le 

premier etage d' egalisation est agence pour produire le 
signal intermediaire sous la forme d'un vecteur Y' de n+p 
echantillons <Y' x , . . . ,y' n+p > obtenu selon la relation : 

25 Y' = <A' H A' )- 1 A' H Y 

ou Y est un vecteur forme de n+W echantillons (y k ) du 
signal re^u, et A' est une ma trice de n+W lignes et n+p 
colonnes ayant une structure de Toeplitz formee a partir 
des coefficients (s q ) dudit polynome en Z" 1 de degre W-p 

30 (R S (Z)). 



+ 
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7. Recepteur selon la revendication 5 ou 6, dans 
lequel le second etage d' egalisation est agence pour 
mettre en oeuvre un algorithme de Viterbi . 

8. Recepteur selon 1'une quelconque des 
5 revendications 5 a 7, dans lequel les moyens de 

repartition des racines utilisent, pour repartir les W 
racines entre les premier et second ensembles, un critere 
de distance au cercle unite s'exprimant par une distance 5 q 

de la forme 5 q = 1 - |a q | si |a q | < 1, et de la forme 

10 8 q = 1 - l/|<* q | si |ct q > 1, pour l<q<W. 
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